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CT. Abstract 

We look for "static" spherically symmetric solutions of Einstein's 
I— li Equations for perfect fluid source with equation of state p = wp. In 

order to include the possibilities of recently popularized dark energy 
^ ■ and phantom energy possibly pervading the spacetime, we put no con- 

, straints on the constant w. We consider all four cases compatible with 

I the standard ansatz for the line element, discussed in previous work. 

■ For each case we derive the equation obeyed by the mass function or its 
, analogs. For these equations, we find all finite-polynomial solutions, 

including possible negative powers. 

For the standard case, we find no significantly new solutions, but 
show that one solution is a static phantom solution, another a black 
I hole-like solution. For the dynamic and/or tachyonic cases we find, 

^ ' among others, dynamic and static tachyonic solutions, a Kantowski- 

■ Sachs (KS) class phantom solution, another KS-class solution for dark 
energy, and a second black hole-like solution. 

The black hole-like solutions feature segregated normal and tachy- 
onic matter, consistent with the assertion of previous work. In the 
first black hole-like solution, tachyonic matter is inside the horizon, in 
the second, outside. 

The static phantom solution, a limit of an old one, is surpris- 
ing at first, since phantom energy is usually associated with super- 
exponential expansion. The KS-phantom solution stands out since its 
"mass function" is a ninth order polynomial. 
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1 Introduction and Motivation 



Exact solutions of Einstein's Field Equations 

= kT^^ (1) 

are, of course, of interest for various purposes. Since the equations are very 
complicated, to find solutions one often makes simplifying assumptions about 
the left-hand-side and/or the right-hand- side. Popular simplifying assump- 
tions about the left-hand-side include staticity and spherical symmetry. As 
is well known, the use of both assumptions together leads to the ansatz [U 
Sect.23.2] 

ds^ = -B{r)dt^ + A{r)dr'^ + r'^d^? (2) 

for the metric. 

Most-often used simplifying assumptions about the right-hand-side of ([T]) 
are that T^i, represents vacuum (i.e. vanishes) or an electromagnetic field 
or a perfect fluid. For example, the vacuum assumption, together with the 
ansatz ([2]) gives uniquely the Schwarzschild metric, the simplest and best- 
known black hole solution. 

The perfect fluid form of T^^,, the stress-energy- momentum tensor, is 

Tf,u = {p + p)Uf,Uu + P9i,u (3) 

where p and p are the energy density and pressure, respectively, as measured 
by an observer moving with the fluid, and is its four- velocity. The use of 
this T^^ together with ansatz ^ describes the interiors of static spherically 
symmetric stars, for example. But the description ^ is not complete: p and 
p should also be specified as functions of particle number density, tempera- 
ture, etc. One further simplifying assumption, justified under most circum- 
stances, is that there is a relation, called an equation of state f{p,p) = 
between p and p. In cosmology, one usually assumes that the equation of 
state is a proportionality, 

p = wp, (4) 

with e.g. w = describing the matter-dominated (or " pressureless dust") 
case, w = 1/3 the radiation-dominated case, w < —1/3 dark energy, and 
w < —1 phantom energy. The latter two concepts have been introduced into 
cosmology in the last decade [21 [3] , after the discovery of the acceleration of 
the expansion of the universe [U E] . 

Now that a good case exists that the universe might be dominated by 
dark energy, even phantom energy, one should look for exact solutions with 
these sources. In particular, static spherically symmetric solutions would 
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be the easiest to find and might be relevant in the contexts of black holes 
or static stars. These solutions can be found starting from the ansatz 
which for "static" perfect fiuid source, (i.e. = u^Sq) leads to the well- 
known Oppenheimer-Volkoff (OV) equation [6] 

where 

F{r) = n pr'^dr, (6) 



and prime denotes r-derivative. F{r) can be recognized as h/Atc times the 
"mass function" defined in the literature. Into the OV equation one must 
put p in terms of p via the equation of state, then p in terms of F', via 
(E]), eventually getting a differential equation for F. After solving for F, the 
metric functions can be found via 

A(r) = (7) 

V ^ r - F{r) ^ ' 

B'{r) npr'^ + 11 

B{r) ~ r - F(r) ~ r' ^ ' 

The solutions can be interpreted as static only for positive A{r) and B{r), 
however. In general, the ansatz admits four classes of solutions, called 
NS (the standard case), TD, ND (corresponding to Kantowski-Sachs P, 
Sect. 15. 6. 5], [7] case) and TS in [U]. The ND and TD solutions are not 
static, hence the quotes on "static" in the title and abstract. For each class, 
one gets a different OV-hke equation. 

The OV equation is valid in case NS. For equation of state (jlj), it becomes 

{w + l)F'{wrF' + F)+ 2w{rF" - 2F'){r - F) = (9) 

where F{r) is written as F for brevity, and we put no constraint on w other 
than that it is a constant. This is a nonlinear equation whose general solution 
is difficult to find. One can attempt a series solution 

oo 

^W = E«n^" (10) 
n=0 

but the recursion expression one gets for a„ involves all of . . . a„_i and it 
seems not possible to even show that (ITUI) converges, let alone find a closed 
expression for a„. 

We can, however, find all of the finite-polynomial solutions of ([9]). This 
we do in the next section. In fact, we find all finite Laurent polynomials. 
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i.e. we consider also negative power^ of r, but find none in case NS. Four of 
the found solutions are valid for particular values of w, and two for general 
w. While none of the solutions is totally original, the procedure shows that 
there are no other finite-polynomial solutions; and in Section [3] we discuss 
properties of the spacetimes. 

In Section m we discuss similar solutions, derived in the appendix, for the 
TD, ND(KS) and TS cases. We also ask if we can find any solutions with 
finite-polynomial A{r). 

2 All finite-polynomial solutions for the mass 
function from the standard OV equation 

In case NS, any power of r less than 3 in F{r) means a diverging density 
at the origin; in particular, a constant term corresponds to a point mass 
there, while negative powers mean diverging mass function, and therefore 
seem unnatural. On the other hand, the meaning of F{r) is different in the 
TD, ND(KS) and TS cases, therefore negative powers are more acceptable. 

The highest and lowest powers of r in F{r) we will call m and m. The 
second-highest, third-highest, second- lowest and third-lowest powers of r in 
F{r) we will call n, p, n and p respectively, when they exist; and A, B, C, A, 
B, C will be the respective coefficients. We will substitute the polynomial 
into the left-hand-side of © and set coefficients of all powers of r equal to 
zero. 

For m > 1, the highest power of r in ([9]) is 2m — 1, but the expression 
for its coefficient will change if m = 1, because of the r-term in (r — F) (and 
similar first power terms in TD, ND(KS) and TS cases). For m = 0, the 
highest power will be n (if it exists), and for m < 0, it will be m. Similar 
considerations apply to the lowest power of r in ([HD as well. Therefore we 
are led to the matrix of cases shown in Table [H 

For cases 1-4, that is, for m > 1, the highest power of r in eq.Q is 2m — 1, 
with coefficient 

(w + l)m(wm + l)y4^ - 2wm(m - 3)y4^ = (11) 

therefore in these cases A is arbitrary and 

7w + 1 

m = — -. 12 

w[l — w) 



^In the rest of this work, we will use "power" also when we really mean "order of the 
power" . It should be clear from the context which meaning is intended. 
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m > 1 


rh = 1 


m = 


m < [eq.([l5D or (^] 


m > 1 [eq.([I2D or (USD] 


1 


2 


3 


4 


m = 1 




5 


6 


7 


m = 






8 


9 


m < 
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Table 1: Matrix of cases for solution of equations and (jBl) by finite 
Laurent polynomials with highest power m and lowest power rh. In the first 
row and last column, an important equation valid for that row/column is 
indicated. 



If this had been an integer, we would have found the order of the poly- 
nomial for arbitrary w. Since it is not, we conclude that in these cases finite 
polynomial solutions exist for certain values of w only. 

Of course, one can also solve for w in terms of m: 



m — 7 ± J {m — 7)2 — 4m 



or write 

= . (14) 

m 

Similarly, for m < (that is, for cases 4,7,9 and 10), we can consider 
lowest power of r in eq.Q and find 

^ = ~7r^~T- ^^^^ 

w[l — w) 

Now, we can start the separate consideration of the cases in Table [U 



Case 1. m > 1, m > 1 



For m > 1, the lowest power of r in eq.dH]) is m, the contribution coming 
from the right part. Its coefficient is 

2w[m{m — 1) — 2m] (16) 

w = is incompatible with eq. (fTTl) . therefore 

m = 3. (17) 

We would like to consider the second- highest power of r in eq.([n]) now. 
The form of that coefficient will depend on the existence of n: 
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Case 1.1. m > 3, m = 3. 

In this case, n exists and the second-highest power of r in eq.Q is m + 
n — 1, with coefficient 

{w + l)[m{wn + l)+n{wm+l)]AB -2w[m{m-3)+n{n-3)]AB = 0. (18) 

giving arbitrary B and, after ehmination of w"^ by using ( |T4l) . the equation 

(4mn + 7m — In — 2m^ — 2n^)w + m — n = (19) 

which not only gives n in terms of m and but also means that w is rational. 

A careful inspection of ( |T3|) shows that there are only three values of m 
giving rational w: 18, 15 and 3. For each m, there are two w values, making 
a total of four subcases of subcase 1.1, since m ^ 3: 

Case 1.1.1. m = IS, n > 1, w = \. 

In this case, solving ffT^ for n gives the values 18 and 27/2, both of 
which are unacceptable; the former because we should have m > n, the 
latter because it is not an integer. Hence, this case fails. 

One can also see this failure using a 'brute force' approach: If one puts 
a general 18**^ order polynomial (in effect, extending the argument down 
to m = 0) into the left-hand-side of ([9]) and sets the coefficients of pow- 
ers of r to zero, starting from the highest (35**^) power, one gets Oig = A, 
an = aiQ = ... = 02 = 0, ai = oq = by the time one arrives at the 17*^^ 
power. But when this polynomial is put afresh into the left-hand-side of (Q, 
one gets instead of zero, so lower powers don't cancel entirely. This is 
not surprising, since there are 36 powers of r in ([9]), but 19 coefficients to be 
found. 

Case 1.1.2. m = 18, n > 1, w = |. 

This time, for n we get 18 and 10, so we should take the latter. We could 
then continue, separating out the third highest power, but the 'brute force' 
approach is more straightforward, especially since it can be executed with 
software. We find that this case also fails. 

Similarly, we find that 

Case 1.1.3. m = lb, n > 1, w = ^ 
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and 



Case 1.1.4- m = 15, n > 1, w = ^ 
fail too, finishing subcase 1.1. 

Case 1.2. m = rh = 3. 

Tliis subcase gives us two solutions, 



Solution 1 
Solution 2 



w 
w 



-1, F(r) 
-I Fir) 



-- Ar\ 



(20) 
(21) 



which finish case 1. 

Case 2. m > 1, m = 1 

The lowest power in Q is now 1 , the vanishing of whose coefficient gives 

4w 



A 



w'^ + 6w + 1 



22 



unless w = —3 ± 2^/2 (For these values, the coefficient cannot vanish at all). 
For more information, we consider the second-highest power in (Q, m + n — 1. 
n exists, but it may or may not be equal to m = 1. This necessitates consid- 
eration of two subcases: 

Case 2.1. m > 1, m = l,n > 1 

In this subcase, eq.(|T8l) is again vahd, therefore the same chain of argu- 
ments can be followed ending with rational w and allowed m values of 18, 
15 and 3. Subcases 2.1.1 - 2.1.4 [fn = 18, 15) are covered by subcases 1.1.1 
- 1.1.4 ill the 'brute force' approach, giving no solutions. This leaves 

Case 2.1.5. m = 3, n > 1, w = —1, 

which fails because the solutions for n are 3 and 0, 



Case 2.1.6. m = 3, n > 1, w = — |, 
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which fails because solutions for n are 3 and 1; this finishes subcase 2.1. 

Case 2.2. m > 1, m = 1 = n 

Putting F{r) = Ar"" + Ar into the OV eqn. ([9]), and using (fT2l) and (l22l). 
we get 

(1 + 3w) (w^ + 6w + 1) = 0. (23) 
As pointed out after eq. (l^ . {w^ + 6w + 1) cannot vanish , so we get 

1 3 
Solution 3 : w = — , F{r) = Ar^ + -r. (24) 

3 2 

This finishes case 2. Solution 3 does not include Solution 2 as a special 
case. 

Case 3. m > 1, m = 

This time, the lowest power in is n — 1, with coefficient 

[{w + l)n - 2wh{h - 3)]AB (25) 

hence 

which again means that w is rational and leads to the same m — w pairs as 
in subcase 1.1. The subcases 3.1 - 3.4 {m = 18,15) are again covered by 
subcases 1.1.1 - 1.1. 4 in the 'brute force' approach with no solutions, leaving 
the m = 3 cases. 

Case 3.5. m = 3, m = 0, w = — 1. 

This leads to fi = 3, therefore to 

Solution 4 : w = -I, F{r) = Ar^ + C (27) 

which includes solution 1 ClS db special case. 

Case 3.6. m = 3, rh = 0, w = — |. 

This gives n = 2, but also n = 1 (see subcase 2.1.6 or eq. ffT^ ). an impos- 
sibility. This concludes case 3. 



8 



Case ^. m > 1, m < 

In this case, eqs. ( fT2l) and ( fT5l) mean that m = rh, hence this case fails. 
Case 5. m = 1 = m 

Very straightforwardly, one can derive 

Aw 

Solution 5 : w arbitrary, F{r) = — r (28) 

+ 6w + 1 

Note that the coefficient A was arbitrary in solutions 1-4, but it is deter- 
mined in terms of w in solution 5. Also, solutions 1-4 required certain values 
of w, while solution 5 is valid for arbitrary w. 

Case 6. m = 1, m = 

It turns out that for one value of w, one can add a constant term to the 
above solution: 

Solution 6 : w = --, Fir) = 5r + B. (29) 

5 

Case 7. m = 1, rh < 

In a sense, this case is a mirror image of case 2. We have (1151) . and 
consideration of highest power of r in ([9]) gives 

A = . (30) 

As in case 2, [tv^ + Qw + 1) cannot vanish. For more information, we con- 
sider the second-lowest power fh+n—1, distinguishing if n is equal to 1 or not. 

Case 7.1. m = 1, rh < 0,h < 1 

The coefficient of r'"+"~^ is given by the same expression as eq. fll8l) with 
m ^ fh, n ^ n, A ^ A and B ^ B. This makes again w rational, but now 
rh must be 18 or 15 or 3, unacceptable because they are positive. 

Case 7.2. m = 1, rh < 0,h = 1 

F{r) consists of two terms, F = Ar + Ar"^ now. The vanishing of the 
coefficient of r"^ in ([9]) reduces upon the substitutions (ITSll and ( |30l) to ( l23ll 
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again, giving the unacceptable (positive) fh value 3. 
Case 8. m = 0, = rh 

This case is trivial: 

Solution 7 : w arbitrary, F(r) = A. (31) 
Case 9. m = 0, m < 

The highest power in ([H]) is n now, with coefficient 2wn{n — 3)B, similar 
to lowest power in case 1. This cannot vanish; hence there is no solution in 
this case. 

Case 10. m,m < 

The same argument as above is valid here for n m, so again there is 
no solution. 

This completes all finite polynomial solutions of equation ([9]). Since So- 
lution 1 is a special case of Solution 4, we will not consider it separately in 
the following section. 

3 Discussion of the solutions found from the 
standard (NS) OV equation 

To finalize the solutions, we calculate the metric functions A{r) and B{r) 
by using ([7]), (IHl), (H]) and (EI). The calculation of B{r) involves an arbitrary 
multiplicative constant at the last stage, the change of which is usually inter- 
preted as a rescaling of t, therefore physically irrelevant. But such rescaling 
cannot change the sign of that constant, so we consider the two choices of 
sign as two separate solutions, unless the requirement of correct signature 
forces a choice upon us. This happens for solutions 3, 4 and 7, whereas for 
solutions 2, 5 and 6 we have consider both signs. The results are shown in 
Table [2l where the well-known solutions are indicated in italics. 

When the metric functions are negative, the spacetime cannot be sup- 
ported by normal perfect fiuid, the source fiuid must be tachyonic. In other 
words, such a spacetime is of type TD in the terminology of [9J. In that 
case, the OV equation, ([5]), is not valid, but still, A{r)-B{r) pairs satisfy the 
same equation of pressure isotropy for cases NS and TD. Therefore negative 
metric functions found from NS-equations represent a valid TD solution, but 
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Sol.No. 


w 


Fir) 


Bir) = -gtt 


y4(r) = fifrr 


Comments 


2 a,b 


1 
3 




±1 


i 


2a, A > : ESU; 

2a, A < : open, static; 

2b : type TD 


3 


1 

3 


Ar^ + |r 




1 


A>0: type TD; 
A<0 : BH-like 


4 


-1 


Ar^ + C 


1 _ C _ ^^2 

r 


i 


/Cd'tt/er (SdS) 


5 a,b 


arbitrary, 
except —1, 
-3 ±2^2 




±(^)^ 




5a: type NS, 

incl. static phantom; 

5b: type TD 


6 a,b 


1 

5 


5r + 5 


±10 

r 


1 

r 


6a: type NS; 
6b: type TD 


7 


arbitrary 


A 


(i-f) 


i 

r 


Schwarzschild 



Table 2: All finite-polynomial solutions of the equation ^ for the mass 
function in the standard (NS) OV case, together with the corresponding 
metric functions. Although we started with the NS OV equation, some of 
the solutions belong to class TD, as defined in [9]. In Solutions 2, 5 and 6, 
the upper signs in B{r) apply to solutions a and lower signs to solutions b. 
The well-known solutions are indicated in italics. 
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not with the equation of state that one has started with. If the NS equation 
of state is the corresponding TD equation of state becomes p = —j^^P- 

Solution 7 is the Schwarzschild solution. It may at first seem surprising 
that there is no restriction on w. But since p vanishes, the value of w does 
not matter. In other words, it corresponds to a situation where all the fluid 
-whatever its equation of state parameter is- has already collapsed to the 
origin. Also, here we do not apply the usual restriction that A must be 
positive. If A is negative, the spacetime will give a naked singularity. 

Solution 2a with positive A is also well-known: it is the Einstein static 
universe, with intimate historical connection to the cosmological constant 
A, equivalent to w = —1. But this universe also contains matter {w = 0), 
whose attraction is precisely balanced by the repulsion of A. So the matter 
density is proportional to A and the net effect is equivalent to a single fluid 
with w = — |. Of course, "in the universe" Ar"^ < 1, so A(r) is positive and 
the signature correct. For negative A, Solution 2a represents an open static 
universe, albeit with negative energy density, and no coordinate restriction. 

Noting that the third well-known solution in Table [5] is Solution 4, the 
Kottler (aka Schwarzschild- de Sitter) solution, the de Sitter part sometimes 
being called anti-de Sitter if A is negative, we turn to the discussion of less 
well-known solutions; changing the order slightly in the interest of presenta- 
tion. 

Solution 2b : 

This solution has correct signature only for Ar"^ > 1, which means A 
must be positive. It is a dynamic spacetime, r being timelike, (it is solution 
TDl of in]) and describes a spacetime that first contracts, then expands in 
angular directions, while distances in the orthogonal spacelike direction stay 
fixed. Even though we found this solution for w = — |, the equation of state 
is actually p = p- The solution can be identified with the n = choice of 
Tolman V [I3], if const=-l is choseiif], with identification — > A. 

Solution 5a : 

Correct signature means positivity of A{r) in this solution, which in turn 
means that the solution is valid except for —3 — 2-\/2 < w < —3 + 2-\/2 (and 
it is of type NS). The cases w < — 3 — 2a/2, for example, w = —6, represent 
static (ultra) phantom solutions. The w = | case is well-known Prob. 

^The KS-like form of the metric is ds'^ = -dr'^ + dp'^ + \ cosh^(%/34r) dft'^. 
^But Tolman chooses const=i?^ and later hterature reports this form (e.g. [T4]). 
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23.10]; the w —>■ oo limit, meaning zero density but nonzero pressure, is the 
metric called SI in [TU]; other valid cases with integer power of r in B[r) are 
w = 1 and w = 3. 

The density is proportional to ^, but this is a mild singularity because 
the mass function goes to zero as r — >■ 0, i.e. there is no mass point at the 
origin. Of course, there is no event horizon, so the singularity is naked. 

The origin is attractive to test particles for w < —3 — 2\/2 and for w > 0; 
repulsive for —3 + 2^/2 < w < 0. The sign of attraction correlates with the 
sign of p + 3p, the so-called "density of active gravitational mass" (e.g. [TT] ) 
for the fluid. The pressure is positive for all w ranges, and since p oc p oc ^, 
the pressure gradient is always towards the origin. So, we cannot understand 
the balance of a fluid element as in terms of p (both forces would accelerate 
the fluid element towards the origin in the ultraphantom case), but in terms 
of p + p, the "density of inertial mass" (e.g. pT]). 

These static ultraphantom solutions constitute a counterexample to the 
impression in the literature (e.g. see [T2j) that everywhere-phantom static 
spherically symmetric solutions cannot exist. 

This solution can be identifled with the n = R oo and B = — 

w+l ' ro 

(or const=rQ ™^^) choice of Tolman V [T3] . 
Solution 5b : 

This is a TD solution (a subcas^ of TD2 of ^) valid for —3 — 2\/2 < 
w < —3 + 2v^, except w = —1. Assuming r is future-directed, this space- 
time expands in the angular directions, and either expands (for w < —1) or 
contracts (for w > —1) in the orthogonal spaoehke iirectio,! An inflnite 
number of tf- values, crowding -1, exist that give integer power of r in B{r). 

The equation of state is p = —j^^P- The density is still proportional 
to ^, but because of the timehke nature of r, F(r) cannot be interpreted as 
the mass function, and therefore we cannot make the same claim as to the 
mildness of the singularity as in solution 5a. 

This solution can be identifled with the almost same subcase of Tolman 
V [13] as solution 5a, except^ const=— rg 

Solution 3 : 



^Which subcase it is depends on the sign of w 

I T, a I uu TTi , vvneie — , , , Yi 



^Metric in KS-likc form: ds^ = -dr^ + ( ™^ dp^ + df^^^ where A - +'^"'+^ 
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For positive A, this solution is also of type TD, contracting in the angular 
directions and expanding in the orthogonal spacelike directioij^ as r — > 0. 

For negative A, both metric functions switch sign at r = rn = 
that the spacetime is static (NS) for r > th and dynamic (TD) for r < r^- 
As far as test particle motion is concerned, this spacetime would be that of 
a black hole; but it must be supported by normal matter in the NS region, 
and tachyonic matter (with p = p) in the TD region. As unreasonable as 
this may seem, it is the only possible perfect fluid interpretation [9]. 

As the origin of r is approached, the density again diverges like ^2, but 
again r is timelike near the origin, and the same (non) conclusion applies to 
the singularity as in Solution 5b. 

For positive A, this solution can be identified with the a = 1, 6 = — 1, 
m = 0, ^ = A (and the trivial i? = 1 or const=l) choice of Tolman VIII [13] . 

Solution 6a : 



This solution is type NS. C must be positive and r < j. Interestingly, 
radially moving free particles oscillate between a minimum radius and j , 
which may be understood in terms of the repulsion of the negative mass point 
at the origin (C = —B and F{r) is the mass function) versus the attraction 
of the fluid, whose "enclosed active gravitational mass" (e.g. [TT]) grows with 
r (here, both p and p + 3p are positive). 

The origin is a naked singularity, and not only due to the negative point 
mass there: The scalar curvature is ^, that is, it diverges without con- 
taining C. But, after all, the scalar curvature does not contain M in the 
Schwarzschild case, either (in fact, it vanishes), r = j is a type of boundary, 
it is a turning point for all radial timelike geodesies. 

This solution can be identified with the n = — |, R ^ —C and B"^ = vq 
choice of Tolman V [13] . 



Solution 6b : 



This TD solution (with p = p/3) can be identified with the n 



_ 1 

2' 

R ^ —C and const =-ro choice^ of Tolman V [13]. There is no coordinate 
restriction for negative C, but r must be larger than j for positive C. 

In the latter case, again r = j is a turning point for timelike radial 
geodesies, but r is timelike, so this spacetime first contracts in the angular 
directions while expanding in the orthogonal spacelike direction, then the 
evolution reverses. 



^The KS-likc form of the metric is ds^ = -dT'^+Acoih^{^/AT)dp'^ + ^ sinh^(V^T) dVt 
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On the other hand, for negative C, the spacetime expands in the angular 
directions while contracting in the orthogonal spacelike directior0, assuming 
r is future-directed. 



4 Discussion of solutions found from the OV- 
like equations in the other cases and an- 
other attempt 

As discussed in the previous section, the TD solutions satisfy the same equa- 
tion of pressure isotropy as the NS solutions, therefore a solution derived 
from the OV(like) equations of one class may in fact belong to the other 
class. Moreover, the same is valid for the ND and TS classes. 

It turns out that for the equation of state p = wp, the OV-like equations 
of the TD and TS cases do not give any solutions not already covered by 
NS and ND cases, except for the special w value — |. This should not be 
taken as an indication that the TD and TS solutions are trivial relabelings; 
for more complicated equations of state, there will be different solutions. 

The proof of the above statement, the application of the procedure of 
Sect. 2 to the ND case to find all ND and TS solutions with finite-polynomial 
F{r) for w 7^ — i, and the TD and TS solutions for w = are given in the 
appendix. 

In this section, we calculate the metric functions A{r) and B{r) for each 
solution from the appendix by using the relevant formulae, and discuss the 
solutions. We also show that no nontrivial solution with finite-polynomial 
A{r) exists. 

4.1 The TD case 

For Solution 8, we get 

which for r < C (only possible if C > 0) gives 
B = -rf^ 

i.e. Solution TD3 of [9]. 



(2r2 + 5Cr - 15C^) + J^-^ ( - IbCHan-' J^-^ 



(33) 



'^KS-likc form of the metric is ds^ = -dr^ + -4p-dp'^ +r^{T) dfl'^ , where ^ = ±W4 - ^ 
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On the other hand, for r > C we find 



B 



(2r2 + 5Cr - 15C^ 



IC-r 



Ci + ln( 



2, ,Vr-C + v^. 



(34) 



the solution called NSl in pj, found in [15] and named Kuch68 I in [13]. It 
describes a spacetime where pure pressure is in static equilibrium with its 
own gravitational attraction. 



4.2 The ND(KS) and TS cases 

The solutions found for the ND(KS) and TS cases, together with their metric 
functions, are shown in Table [3] (Solution 9 does not appear because it is a 
special case of Solution 11). As in Sect. 3, the sign of B{r) is arbitrary, unless 
forced by the signature requirement. 

The Schwarzschild and Kottler (SdS) solutions, which appeared in Table 
[21 are found in this table as well, because they cannot really be classified 
in this scheme. Our classification is based upon the nature and direction 
of motion of the fiuid, but for these solutions, the stress-energy-momentum 
tensor is independent of the fiuid four-velocity: The u^u^ term in T^i, is 
multiplied by p + p; and p + p = for the Kottler solution, p = p = for 
Schwarzschild. Hence, these solutions satisfy the equations for all four cases. 

The other solutions in the table are less well-known: 



Solution 10 : 



For positive A, this solution is type ND (KS), representing a dynamic 
spacetime filled with a phantom perfect fiuid. Assuming r is future-directed, 
the spacetime expands in the angular directions; in the perpendicular space- 
like direction, it first contracts, reaches a minimum, then expand^. It is 
singular at both ends of the evolution, that is, at r = and as r — oo, the 
first singularity being in the finite past, the second in the infinite future. Of 
course, these attributes switch if r is past-directed. 

For negative A, Solution 10, like Solution 3, represents a black hole space- 
time, as far as test particle motion is concerned; but it must be supported by 
two different fluids on the two sides of the horizon: tachyonic fluid in the out- 
side, static region and normal fluid in the dynamic region inside/in the future. 



^The KS form of the metric is ds^ = -rfr^ + l+^^ll^dp^ + r'^{T)dn^, where 
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Sol.No. 


w 


F(r) 


B{r) = -gu 


v4(r) = 


Comments 


10 


-3 




L-\-0/il 

^6 


Q 
O 

l+8Ar« 


A > 0: type ND(KS 
phantom-filled dy- 
namic universe; 
A < : BH-like 


11 


-1 


Ar^ + B 


1 - ^ - ylr2 

r 


J 

l-^-Ar^ 


Kdttler (SdS) 


12 a,b 


arbitrary, 
except 
— 1 and 1 


Sw^— 2ii;— 1 


T 


(w-l)(3w+l) 
(«;+l)^ 


— 1 < w < 1: type 
TS; otherwise: 
type ND(KS), inch 
DE, inch phantom 


13 a,b 


1 

3 


-|r + i? 


±10 

r 


4r-3B 


13a: type TS; 
13b: type ND(KS) 


14 


arbitrary 


A 




1 

1-4 


Schwarzschild 


15 a,b 


1 


c_ 

r 


±1 





15a: type TS; 
15b: type ND(KS) 


16 


1 

2 


c 




1 

1-^ 

r 


Schwarzschild 


17 


1 
2 


4 

3^ 


_/Jl\-4 


-3 


type ND(KS) 



Table 3: All finite-polynomial solutions for F{r) in the ND(KS) and TS 
cases as defined in [9]; together with the corresponding metric functions. In 
Solutions 12, 13 and 15, the upper signs in B{r) apply to solutions a and 
lower signs to solutions b. The well-known solutions are indicated in italics. 
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Solutions 12a,h : 

Solution 12a is a TS solution, valid for — | < w < 1. For positive w, 
radially incoming test particles are reflected near the origin back to infinity, 
whereas for negative the origin constitutes a potential well from which 
they cannot escape. 

Solution 12h, however, is valid for w values other than — ^ < w < 1, and 
is identical to the Ci = special case of Solution ND2 of [9J. If r is future- 
directed, it expands in the angular directions, and it expands or contracts in 
the perpendicular spacelike direction, if the sign of ^^^^ is negative or pos- 
itive, respectivel}{§. Note that this means expansion for non-phantom dark 
energy (—1 < w < — |) and "radial" contraction for phantom energy. 

Solutions 13a,h : 

Solution 13a is a TS solution, where we must have 4r > 35, i.e. we have 
a restriction on r if i? is positive. Either way, the equation of motion for test 
particles shows that tachyonic w = | fluid is repulsive, consistent with the 
solution 12a. 

Solution 13h is an ND (KS) solution, where 4r < 3i?. It represents a 
radiation-filled universe that expands and recoUapses in angular directions 
while contracting and reexpanding in the perpendicular spacelike directiorF^: 
first found in [TB] . 

Solutions 15a,h : 

Solution 15a is solution TSl of [9J, where we must have > C. For 
positive C, tq = VC is a turning point for radial geodesies; for negative C, 
there are no such turning points. 

Solution 15b is solution NDl of [9], apparently first found in [17], describ- 
ing a finite-lifetime universe containing stiff matter, expanding and recoUaps- 
ing in the angular direction^. 

Solution 17 : 



^The KS form of the metric is ds"^ — —dr^ + j '"^^ dp^ + j^dil^, where 

< _ (t»-l)(3t»+l) 
^ - + • 

10 The KS form of the metric is ds^ = -dT^ + :;^^dp^+r^{T) dn^, where ^ = ±y^f - |. 

The arbitrary ro must be chosen as jB for agreement with [16J, p. 1684. 
"The KS form of the metric is ds^ = -dr^ + dp^ + (C - t^) dn^. 
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This solution is the Ci = 0,A = —3 special case of solution ND2 of 
[9], describing a spacetime containing pressure, but no density (because it 
is an ND (KS) solution found from the TS equations, its equation of state 
is not p = — |); expanding in ang ular directions while contracting in the 
perpendicular spacehke directioio, if r is taken to be future-directed. 

4.3 Finite-polynomial A{r)? 

Another possible way to look for solutions is to work in terms of A{r) rather 
than F{r) by using equation ([7]). This leads to an equation with terms 
second to fourth order in A{r) and/or its derivatives. In trying to find a 
finite-polynomial solution for A{r), if the highest power of r in A{r) is m, the 
highest power of r in the equation is 4m; but it is multiplied by A^^w + lY in 
cases NS and TD, and —A'^{w + 1)^ in cases ND and TS. Setting the trivial 
w = —1 case aside, therefore, the highest possible value for m is zero. A 
similar argument shows that the lowest power in the A{r) polynomial must 
be zero or higher. Hence the only finite polynomial A{r) can be for equation 
of state (jll) is a constant. 

5 Summary and Conclusions 

We have considered spherically symmetric perfect fluid solutions in Gen- 
eral Relativity and found all finite-polynomial solutions -including negative 
powers- of the equation satisfied by the so-called "mass function" and its 
mathematical analogs for the equation of state p = wp; and discussed the 
associated spacetimes. 

The equation for the mass function follows from the Oppenheimer-Volkoff 
(OV) equation in the standard case where the fluid is static and normal (i.e. 
non-tachyonic, m'^m^ = — 1). However, the metric ansatz used in that analysis 
can also accomodate cases where the spacetime is dynamic in a certain way, 
or the fluid is tachyonic; as discussed in [S]. In these other cases analogous, 
but different functions exist, satisfying their own equations. 

The solutions we found for the standard case, NS, are mathematically not 
very original; they are either some limiting cases of solutions found long ago 
by Tolman |13] or simple modifications thereof. Some aspects of the physical 
nature of these solutions can be seen in new light however, considering the 
classification in P] and newly cosmologically relevant concepts of dark energy 

12 The KS form of the metric is ds^ = -dr^ + dp^ + ^ dVL^ . 
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and phantom energy. The solutions (Table [2]) include dynamic spacetimes 
supported by tachyonic fluids {2b, 5b, 3 with A > 0, 6b) and a static space- 
time containing a w = — | fluid around a negative point mass {6a). The 
TD case gives two extra solutions, one describing a spacetime where pure 
pressure is in static equilibrium with its own gravitational attraction. 

Some interesting solutions are also found from the ND(KS) and TS cases 
(Table ED: There are static solutions supported by tachyonic ?l\xids {12a, 13a, 
15a), the flrst two presumably original. Some solutions {10 for positive A, 
12b, 13b, 15b, 17) are of the Kantowski-Sachs (KS) class: Solutions 13b, 15b 
and 1 7 describe dynamic KS-universes containing radiation, stiff matter and 
pure pressure, respectively. 

We would like to particularly point out the following solutions: 

• Solution 5a for w < —3 — 2\/2 represents, perhaps unexpectedly, a 
family of static "ultraphantom" solutions. 

• Solution 3 for negative A is a black hole-like spacetime, which must be 
supported by normal matter outside the horizon and tachyonic fluid on 
the inside. 

• Solution 10 for positive A is a phantom KS solution, probably new. 

• Solution 12b can also be valid for dark energy, including phantom, 
exhibiting anisotropic expansion for non-phantom dark energy. 

• Solution 10 for negative A is similar to Solution 3, a black hole-like 
spacetime, supported by segregated normal and tachyonic matter, ex- 
cept in this solution, the tachyonic fluid is outside and normal fluid is 
inside. It was concluded in [9] that black holes supported by perfect 
fluids cannot be "simple". 

There are no other solutions where F{r) is a flnite polynomial of r for the 
assumed equation of state. One can also express the problem(s) in terms of 
A{r), and then try to flnd flnite polynomial solutions. The only such solution 
is yl(r)=constant. 
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A All finite-polynomial solutions for F(r) from 
the OV-like equations in the TD, ND(KS) 
and TS cases 

A.l The TD case 

The TD OV equation [9J is 

, (/tpr^ + Ftp) ^ . , . 

^= 2rir-Frn) + ^''^ 

where 

FToir) = -kJ{p + 2p)r^dr, (36) 



and the metric functions are found by 



r 



A = 77- (37) 



B' Kpr^ + 1 1 



(3^ 



B r — Ftd f 

The substitution p = — (p + 2p) brings the TD OV equation ( l35l) into 
the same form as the regular one in terms of p and p. When expressed 
in terms of Ftd, with equation of state (jlj), we get equation ©, but with 
the replacement w — > —j^^- Since this is another constant equation of 
state parameter, we will not get any finite-polynomial solutions that are not 
already in Tabled unless 1 + 2w = 0. In that case, Ftd becomes a constant. 

Solutions: w = FTD{r)=C. (39) 

A. 2 The ND case 

The ND OV equation [9] is 

, 3FjvD - 4r + Kpr^ 

P= — IT-, ^ — ^ — P + P 40 

2r(r - F^d) 

where 

FNoi't') = J pr'^dr (41) 

and the metric functions can be found by 

r 

A = TT- (42) 

t-Fnd ^ ^ 



B' _ 1- Kpr^ 1 
B r - Fnd r 



(43) 
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In this case, Fj^D^r) obeys 
{l + w)F'j^j,{3wFND-^wr-rF'j,j,) + 2w{rF'j:,j,-2F'j^j,){FND-r) =0 (44) 

To find all finite-polynomial solutions of this equation (dropping label 
ND), we follow the same procedure as in Section 2, and for similar reasons, 
we use the same matrix of cases, shown in Table [H 

For cases 1-4, the consideration of highest power of r in eq.(jllj) gives 

{w-l){3w + m) = (45) 
Similarly, for cases 4,7,9 and 10, consideration of lowest power leads to 

{w -l){3w + m) = 0. (46) 

Case 1. m, m > 1 

Lowest power in eq. (H^ gives 

m = l-2w. (47) 
Two subcases, w = 1 and m = —3w follow from eq. fl45p . 

Case 1.1. m, rh> l;w = 1 

w = 1 gives m = —1 by eq. fHTl) . which contradicts the definition of case 1. 
Case 1.2. m, m > l;m = —3w 

m = —3w gives 

m 

We may have m = rh or m > m: 
Case 1.2.1. m = —3w = m > 1 
gives 

Solution 9 : w = -I, F{r) = Ar^ (49) 
Case 1.2.2. m = —3w > rh > 1 

When n exists and is larger than 1, as in this case, consideration of the 
second-highest power in eq. (H^ gives 

3(m + n)w'^ + [2(m^ + n^) — 2nm — 3(m + n)]w — 2nm = 0, (50) 



im-l] 



(4J 
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which together with m = —3w yields 

m = 2n + 3. (51) 
But, this equation, together with eq. lHHl) gives 



4 

impossible because it violates n > m for m > 1 
Case 2. m > 1, rh = 1 

Lowest power in eq. fH^ gives 



3m — 9 , , 

n=^—, (52) 



A = —^-^ . (53) 

3^2 - 2u; - 1 ^ ' 

unless (3ty2 — 2w — 1) vanishes, that is, w = 1 or w = — |. Highest power 
gives eq. fl45p . so this means m = —3w. n exists, second-highest power gives 
eq-dSni), «/n > 1. 

Case 2.1. m = —3w > 1, m = 1, n > 1 

Again, eq. fl3U|) . together with m = —3w gives eq. flFIj) . Then the third- 
highest power in (jH]) is m + p — 1 = 2r;, + p + 2, as opposed to 2n — 1. We 
have to also distinguish if p > 1 or not. 

Case 2.1.1. m = —3w > 1, m = 1, n > l,p > 1 

Setting the relevant coefficient equal to zero gives 

{m - p){2p + 3 - m) = 0, (54) 
whose solutions, p = m and p = = n are both unacceptable. 

Case 2.1.2. m > 1, m = 1, m = —3w, n > l,p = 1 

Setting the coefficient of the third-highest power in (jUj) equal to zero 
gives 

^ _ ■n{2n + 3) _ 6w{w + 1) 

~ 3(n2 - 1) ~ 9w;2 + 18w + 5 ^ ' 



which contradicts eq.(]53 
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Case 2.2. m > 1, m = 1, m = —?>w, n = 1 

Putting F{r) = Ar"^ + Ar into eq.f H^ . and using (l53|) . followed by 
m = —3w, we get 

(m-9)(m- 1) = 0. (56) 

This gives 

9 

Solution 10 : w = -3, F(r) = Ar^ + -r. (57) 

8 

Case 3. m > 1, m = 

From consideration of the lowest power in ( H4l) . n — 1, we get 

3{w + l) = 2{n-3), (58) 
n exists and according to eq. (H5ll we have two subcases. 

Case 3.1. m> 1, rh = 0, w = 1 

w = 1 gives n = 6 according to fl3Sl) . Now, either n > 6, or n = m = 0. 

Case 3.1.1. m>l, m = 0, w = l, 'ri>6 

In this subcase, eq.( l50|) is applicable, but for w = 1, eq. (!50|) gives 2{m—nY 
on the left-hand-side, which cannot vanish. 



Case 3.1.2. m > 1, m = = n, w = 1, 

When F{r) = Ar^ + A is put into eq. (1441) with w = 1, we get 

72AAr^ - SiAr^ = 0, (59) 

impossible to satisfy. 
Case 3.2. m = -3w > 1, m = 



m = —3w, together with (pa)) gives 



3 — m , , 

n=^— + 3 (60) 



This makes {m,n) pairs (3,3), (5,2), (7,1) possible. 



24 



Case 3.2.1. m = —3w = 3 = n, m = 



Solution 11 : w = -1, F{r) = Ar^ + B. (61) 
This includes Solution 9 as a special case. 

Case 3.2.2. m = -3w = 5, n = 2, m = 

This subcase gives no solution, which can be seen either by the 'brute 
force' method (Sect. 2, case 1.1.1) or by considering the second-highest power 
in eq.( H^ . which gives n = 1 or ri = 5 (fails since we must have 5 > n > 2). 

Case 3.2.3. m = -3w = 7, n = 1, m = 

Again, the 'brute force' method shows that there is no solution. 

Case ^. m > 1, m < 



Since both eq. fHol) and eq. (H6l) apply, we must have 'u; = 1. n exists. If 
n > 1, the argument in case 3.1.1 is applicable, so we must consider n < 1. 
A similar argument for second- lowest power in eq.( H^ gives n > 1, so we 
have F{r) = Ar^ + Br + Ar^. Now, consideration of the coefficients of the 
second-highest and second-lowest powers gives 

m(m + l) m(m + l) 

7 T\2" = B = 7^ -T^ (62) 

(m — ly [m — 1)^ 

with two solutions 

_ m + 1 

m = m, m = (63) 

3m — 1 

neither of which is negative. 
Case 5. m = 1 = m 
Easily: 

Solution 12 : w arbitrary. Fir) = — r. (64) 

As in case ND; this solution is valid for arbitrary w, with A is determined 
in terms of w; whereas in solutions 9-11 A was arbitrary, but w specified. 
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Case 6. m = 1, m = 

As in case ND, one can add a constant term to the above solution for one 
value of w: 

Solution 13 : w = F{r) = — r + B. (65) 

3 3 

Case 7. m = 1, m < 

We have eq. fHUl) . and consideration of highest power in eq. fjl^ gives 



A = — : 66 

3^2 - 2u; - 1 ' ^ ' 

as in case 2, w cannot be 1 or — i, which means rh = —3w. 

Case 7.1. m = 1, m = —3w < 0, fi < 1 

The coefficient of the second-lowest power, r™"^"~^ is given by the same 
expression as in case 1.2,2, with m^rh, n^h, A^A and B B; giving 
us 

m = 2n + 3, (67) 

which means that m < n < —3. We now consider the third-lowest power, 
whose value and coefficient depend on the value of p. 

Case 7.1.1. m = 1, rh = —3w < h < —3, p < 1 

The third-lowest power could be m + p — 1 or 2n — 1. Because of eg. (1671) . 
J9 = — 3 is critical: 

Case 7.1.1.1. m = 1, rh = —3w < n < —3, p < —3 

Power m+p—1 is the third-lowest. Setting its coefficient equal to zero, us- 
ing eg. (1671) and rh = —3w, we get the two solutions p = h and p = 2fi+3 = rh, 
which are both unacceptable. 

Case 7.1.1.2. m = 1, rh = —3w < h < —3, p = —3 

In this subcase, the coefficients A, B and C mix (unlike the previous 
subcase), so we resort to another 'brute force' approach: We put F(r) = 
Ar + B + Cr-^ + Dr'^ + Cr'^ + Br^ + Ar'^^+^ into eq.dSD, use w = 
and sequentially set the coefficients of • • • to zero, for each calculation 
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using the result of the previous one, as well. We get A = g^^2'^"^^-'_^3) (compat- 
ible with eq. (l66i) ). B = 0, C = 0, D = 0, C = in sequence; this last result 
tells us that this subcase fails. 

Case 7.1.1.3. m = 1, m = — 3w <n < —3, p > —3 

Power 2n — 1 is the third-lowest. Setting its coefficient equal to zero, 
using eq. (j6711 and m = —3w, we get h = m, which is unacceptable. 

Case 7.1.2. m = 1, in = —"iw < h < —3, p = 1 

Power 2n — 1 is the third-lowest. Setting its coefficient equal to zero, 
using eq.( l67I) and rh = —3w, we get h = — | or n = — 1 or n = 3 which are 
not acceptable. 

Case 7.2. m = 1, m < 0, n = 1 

F{r) consists of two terms, F = Ar + Ar"^ now. With this F, eq. (H4l) 
reduces, using eq.(l671) and rh = —3w to 

_ 2im^(m-9)(m-l) - 

9(m + 3) ' ^ ^ 

which cannot vanish. 

Case 8. m = 0, = m 

This case is trivial: 

Solution 14 : w arbitrary, F(r) = A. (69) 
Case 9. m = 0, m < 



We have eq. fH6|) . which leads to the subcases w = 1 and m = —3w. The 
highest power in eq. (H^ is n, and consideration of its coefficient gives 



n = l-2w. (70) 

Case 9.1. m = 0, m < 0, w = 1 

Eq. (!70l) gives n = —1, which means that m,h < 1, so the second-lowest 
power in eq. (1441) gives eq. (!50l) . with m —>■ m, n —>■ h. However, for w = 1, 
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the left-hand-side of that equation cannot vanish. 

Case 9.2. m = 0, —3w = rh < 

Since Th,h < 1, eq.( l67|) is vahd. 

Case 9.2.1. m = 0, n = rh = —3w < 

This means that F consists of two terms. We get w = —1 and n = 3, 
unacceptable for this case. 

Case 9.2.2. m = 0, n > m = -3w < 

This means that a p < 1 exists. For p ^ —3, the subcases 7.1.1.1 and 
7.1.1.3 carry over exactly; for p = —3, a calculation analogous to that of 
subcase 7.1.1.2 shows the nonexistence of solutions. 

Case 10. m,fh <{] 

Again, we have eq. (l46l) . therefore subcases w = 1 and rh = —3w. The 
highest power in eq. (1441) is m, consideration of its coefficient giving 

m = l-2w. (71) 

Case 10.1. m,rh <0, w = 1 

Case 10.1.1. m = rh < 0, w = 1 

Solution 15 : w = 1, F{r) = C/r. (72) 
Case 10.1.2. rh < m < 0, w = 1 

n, h exist; subcase 9.1 carries over exactly. 
Case 10.2. m,fh = —3w < 0, 
Case 10.2.1. m = rh = —3w < 

Gives w = —1, m = 3, unacceptable for this case. 
Case 10.2.2. rh = -3w < m < 
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Since m,h < 1, eq. flU7|) is valid. 

Case 10.2.2.1. m = -3w = n < m < 

This means that F consists of two terms. We get w = 5, m = —9 and 
rfi = —15, looking acceptable, but eq. fH^ can still not be satisfied. 

Case 10.2.2.2. fn = -3w < n < m < 

Subcase 5.^.^ applies. 

This completes all finite polynomial solutions of equation (jH]) . 

A. 3 The TS case 

The TS OV equation ^ is 

P+^P= 2r(r - Frs) + ^^'^ 

where 

Frsii") = —1^ j pr'^dr, (74) 

and the metric functions are found by 

r 

A = -T^ (75) 



B' _ 1 + k(p + 2p)r2 1 
B r — Fts r 



(76) 



Again, for Fts, with equation of state (jlj) we get equation (1441) . but with 
the replacement w ~ i+2w ^ leading to the only potentially new solutions 
for 1 + 2w = 0. Then, we get two solutions: 



Solution 16 : w = -i, Frsir) = C, (77) 

4 
3 



Solution 17 : w = -|, Frsir) = -r. (78) 
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